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' We test the vahdhy of the QCD sum rules apphed to the meson ^"''(4430), by considering a 

I diquark-antidiquark type of current with J^ = 0~ and with J^ — 1~ . We find that, with the 

studied currents, it is possible to find an acceptable Borel window. In such a Borel window we have 
simultaneously a good OPE convergence and a pole contribution which is bigger than the continuum 
contribution. We get mz = (4.52 ± 0.09) GeV and mz = (4.84 ± 0.14) GeV for the currents with 
J^ = 0~ and J^ = 1~ respectively. We conclude that the QCD sum rules results favors J^ = 0~ 
quantum numbers for the Z+(4430) meson. 
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During the past years, a series of exotic charmonium like mesons, called X, Y and Z, have been 
discovered in B mesons decays. Among them, the charged resonance state Z+(4430), observed by Belle 
Collaboration [l| in the Z+ ip'-n:^ decay mode, is the most intriguing one since it can not be described 
\ as ordinary cc meson. 

04 ■ The nature of the Z+(4430) meson is cornpletely open and there are already many theoretical inter- 

rn ; pretations about its structure [1 i, i, [E i, 0, i, S [M [Hi [ll E Q E H O ■ However, an 

. intriguing possibiUty is the interpretation as tetraquark or molecular state. In ref. [5|, the closeness of 

O ; the Z+(4430) mass to the threshold of i:)*+(2010)l)i(2420) lead the authors to consider the Z+(4430) 

00 , as a D*Di molecule. This hypothesis was tested in ref. Q by using the QCD sum rules approach, with 

' a good agreement with the experimental data. The interpretation of Z+(4430) as tetraquark state was 

done in refs. [1,H[3. 

Since Z+(4430) was observed in the ip'n^ channel, it is an isovector state with positive G-parity: 
= 1+. However, nothing is known about its spin and parity quantum numbers. For a D* Di molecular 
state in s-wave, the allowed are 0^, 1^ or 2~, although the 2~ assignment is probably suppressed 
in the B — > Z{AA'iQ)K decay, by the small phase space. In this work we use QCD sum rules (QCDSR) 
[20I [2ll . [2^ , to study the two-point function of the state Z+ (4430) considered as a tetraquark state with 
= 0- and JP = 1-. 

In previous calculations, the QCDSR approach was used to study the X(3872) by using a diquark- 
antidiquark current 23], the Z+(4430) meson, by using a D*Di molecular current 9] and the Y mesons 
[2^ by using molecular and diquark-antidiquark type of currents. In all cases a very good agreement 
with the experimental mass was obtained. 

Let us consider first the Z+(4430) by using a diquark-antidiquark current with = 0^ and positive 
G parity. One can invoke simple arguments using constituent quark model to show why the tetraquark 
with the suggested quantum number could be stable. In the constituent quark model, a multiquark exotic 
is expected to have some scalar diquark component in the color anti-triplet configuration, as this is the 
most attractive quark-quark channel. However, when a tetraquark has = 0^ quantum number, it 
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would energetically be more favorable to decay in s-wave into two pseudo-scalar mesons. In terms of 
the spin spin interaction, one can say that the attraction in the quark- antiquark configuration in the 
two pseudo-scalar mesons is phenomenologically more than a factor 3 larger than that in the two scalar 
quark-quark channel in the tetraquark [25|. However, when the tetraquark configuration has = 0~ 
quantum number, at least one of the diquark could be in the attractive channel, while the remaining 
diquark is in the pseudo scalar channel. On the other hand, it can not decay into final states containing 
a pseudo-scalar meson in s-wave; hence the tetraquark could be quasi-stable. To test such configuration 
in a non-perturbative way, we are implementing the QCD sum rule method. 
A possible current describing such state is given by: 

3 = "-^^^[{ulC^,c,){daCcl) - {ulCc,){dai,Ccl)] , (1) 

where the index T means matriz transposition, a, 6, ... are color indices and C is the charge conjugation 
matriz. 

The QCD sum rules for the meson mass are constructed from the two-point correlation function: 

Ii{q) =ijd^x e^'^--(0|r[j(x)jt(0)]|0). (2) 

Phenomenologically, the correlator can be expressed as a dispersion integral 

ds ^ ^ + . . . , (3) 

where pP^'^"{s) is the spectral density and the dots represent subtraction terms. The spectral density is 
described, as usual, as a single sharp pole representing the lowest resonance plus a smooth continuum 
representing higher mass states: 

ff^^^is) - A2^(s-m|)-fp™"*(s), (4) 

where A is proportional to the meson decay constant, fz, which parametrizes the coupling of the current 
to the meson Z~^: 

{0\j\Z+) - fzm% - A. (5) 

It is important to notice that there is no one to one correspondence between the current and the state, 
since the current in Eq. ([T]) can be rewritten in terms of sum a over molecular type currents, by the use 
of the Fierz transformation. However, the parameter A, appearing in Eq. ([S]), gives a measure of the 
strength of the coupling between the current and the state. 

We follow the prescription that the continuum contribution to the spectral density, p'^°"*(s) in Eq. (j4|), 
vanishes bellow a certain continuum threshold sq. Above this threshold, it is given by the result obtained 
with the OPE [26|]: 

*(s) = pO^^(s)e(s-so), (6) 

On the OPE side, we work at leading order in as and consider the contributions of condensates 
up to dimension eight. To keep the charm quark mass finite, we use the momentum-space expression 
for the charm quark propagator. The light quark part of the correlation function is calculated in the 
coordinate-space. Then, the resulting light-quark part is Fourier transformed to the momentum space in 
D dimensions and it is dimensionally regularized at D = 4. The correlation function in the OPE side 
can be written as: 

/■oo OPE( \ 

n°^^((z')=/ ds^ ^+n""<'?''>(q2) , (7) 

J4in2 S — q 

where p'~'^^{s) is given by the imaginary part of the correlation function: Trp^^^{s) — Im[n'^-''^(s)]. 

After equating the two representations of the correlation function, assuming quark-hadron duality, 
making a Borel transform to both sides, and transferring the continuum contribution to the OPE side, 
the sum rule for the pseudoscalar meson up to dimension-eight condensates, is given by: 

J 4mJ 
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The integration Umits are given by amin = (1 ^ -^/l — Am'^/s)/2, amax = (1 + — 4m2/s)/2 and 
/3mjn = am^ / (sa - TO^ ) . 

One should note that a evaluation of the higher dimension condensate contributions is technically 
difficult and non-trivial, which cannot be obtained by a simple routine iteration of the quark propagator 
in an external field. Violation of the factorization hypothesis become increasingly important in higher 
dimensions and so the results become increasingly model dependent, as more condensates will have to be 
introduced if factorization is not valid [l^ • 

Similarly to the results in ref. [o^ , the current in Eq. ([T]) does not get contribution from the quark and 
mixed condensates. This is very different from the OPE behavior obtained for the diquark-antidiquark 
current used for the X(3872) and y(4660) mesons in refs. [23l. [2^. but very similar to the OPE behavior 
obtained for the axial double-charmed meson Tec, also described by a diquark-antidiquark current f2^. 

In the numerical analysis, the input values are taken as 22, 29]: mc{mc) — (1.23 ± 0.05) GeV, {qq) = 
-(0.23 ± 0.03)3 Qgv^ (qga.Gq) = ml{qq) with = 0.8 GeV^ (g^G^) = 0.88 GeV'^. 

We evaluate the sum rule in the Borel range 2.2 < < 3.5 GeV^. To determine the allowed Borel 
window, we analyse the OPE convergence and the pole contribution: the minimum value of the Borel 
mass is fixed by considering the convergence of the OPE, and the maximum value of the Borel mass is 
determined by imposing that the pole contribution must be bigger than the continuum contribution. To 
fix the continuum threshold range we extract the mass from the sum rule, for a given sq, and accept such 
value of So if the obtained mass is around 0.5 GeV smaller than ^/sQ. However, in this case, to be able 
to compare our results with the results obtained by using a molecular type current in ref. j9(], we use the 
same continuum range as in ref. [9]: 4.8 < ^/sq < 5.0 GeV. 

From Fig. [1] we see that we obtain a quite good OPE convergence for M'^ > 2.3 GeV^. Therefore, 
we fix the lower value of in the Borel window as M^j„ = 2.3 GeV^. This figure also shows that 
the dimension-eight condensate contribution is very small as compared with the four-quark condensate 
contribution. 

The comparison between pole and continuum contributions for = 4.9 GeV is shown in Fig. [21 from 
where we see that the pole contribution is bigger than the continuum for IVP < 3.1 GeV^. The same 
analysis for the other values of the continuum threshold gives < 2.9 GeV^ for — 4.8 GeV and 
M2 < 3.3 GeV^ for ^ = 5.0 GeV. 

To extract the mass mz we take the derivative of Eq. ([8]) with respect to and divide the result 

by Eq. In Fig. [H we show the Z+ meson mass, for different values of yso, in the relevant sum rule 
window, with the upper and lower validity limits indicated. From this figure we see that the results are 
very stable as a function of M^. 

To check the dependence of our results with the value of the charm quark mass, we fix ^/s^ — 4.9 GeV 
and vary the charm quark mass in the range rric — (1.23 ± 0.05) GeV. Using 2.5 < < 3.1 GeV^ we 
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FIG. 1: The OPE convergence in the region 2.2 < < 3.5 GeV^ for ^ = 4.9 GeV. Perturbative contribution 
(dotted hne), {g'^G^) contribution (dashed hue), (gg)^ contribution (Ion-dashed Une), {qga.Gq) {qq) (dot-dashed 
line) and the total contribution (solid line). 
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FIG. 2: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution for -ySo ~ 4.9 GeV. 



get: mz = (4.51 ± 0.06) GeV. Including the uncetainty due to the value of the continuum threshold and 
the value of the Borel parameter we arrive at 

mz^^^^ - (4.52 ±0.09) GeV, (12) 

which is a little bigger than the experimental value 1], but still consistent with it, considering the 
uncertanties. Comparing our result with the result obtained in ref. Q: mwDi = (4.40 ± 0.10) GeV, 
where the Z+(4430) was considered by using a D*Di molecular current with = 0^, we see that the 
result in ref. [9| is in a better agreement with the experimental value. However, as mentioned above, 
since there is no one to one corresponde between the structure of the current and the state, we can not 
use this result to conclude that the ^+(4430) is better explained as a molecular state than as a diquark- 
antidiquark state. To get a measure of the coupling between the state and the current, we use Eq. ([8]) to 
evaluate the parameter A, defined in Eq.©. We get: 



Xz^^_^ = (3.75 ±0.48) X 10"^ GeV^ 



(13) 
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FIG. 3: The Z+ with = 0" meson mass as a function of the sum rule parameter (M^) for different values of 
^/sQ■. ^/so = 4.8 GeV dashed line,y'so — 4.9 GeV solid line and ^/so — 5.0 GeV dot-dashed line. The crosses 
indicate the region allowed for the sum rules. 



while for the current used in ref. Q we get: 

Ao.Di = (5.66 ± L26) x 10"^ GeV^ 



(14) 



Therefore, it is possible to conclude that the physical particle with — 0^ and quark content ccud 
couples with a larger strength with the molecular D*Di type current than with the current in Eq.(IT]). 

We now consider the Z+(4430) by using a diquark-antidiquark current with — 1~ and positive G 
parity. The lowest-dimension interpolating operator describing such current is given by: 



^abc^dec r/ T 



V2 



[("a C'75C6)(dd7^75CCe ) -|- (u„ Cj5liJ.Cb)iddj5Cc^ )] 



The two-point correlation function is now given by: 
from where we get 

W^{q) = -iq^n{q% 

and, therefore, we can write a sum rule for n(g'^) as before. The spectral density is now given by 

C^max 1 — a 
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(19) 



Although with this current we stiU do not get contribution from the quark condensate, we do get 
contribution fr'om the mixed condensate. As can be seen by Fig. 31 the mixed condensate contribution is 
of the same order as the four-quark condensate contribution, but with opposite signal. The contribution 
of the dimension-eight condensate is now of the same order as the four-quark condensate contribution, 
for small values of . 
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FIG. 4: The OPE convergence for the sum rule for with = 1~ , using — 5.3 GeV. The dotted, dashed, 
long-dashed, dot-dashed, solid with dots and solid lines give, respectively, the perturbative, gluon condensate, 
mixed condensate, four-quark condensate, dimension- eight condensate and total contributions. 



In this case we find that the continuum threshold is in the range y/sg = (5.3 ± 0.1) GeV and, from 
Fig. m we see that there is a good OPE convergence for > 3.9 GeV^. 

The upper limits for for each value of ^/sq are given in Table I, from where we see that the Borel 
window in this case has higher values of the Borel parameter, as compared with the case for Z+ with 

JP = 0-. 

Table I: Upper limits in the Borel window for Z'^ with — 1~ . 



( GeV) 


Af^^JGeV^) 


5.2 
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5.3 


4.7 


5.4 


5.0 



In the case of Z+ with — 1 we get a worse Borel stability than for the Z+ with = 



the 



allowed sum rule window, as a function of M^, as can be seen by Fig. [51 We also observe that the results 
are, in this case, more sensitive to the values of me- 

Using the Borel window, for each value of sqj to evaluate the mass, and then varying the value of the 
continuum threshold in the range 5.2 < ^/Jq < 5.4 GeV, we get "^z^^.j = (4.80 ± 0.08) GeV. 

Because of the complex spectrum of the exotic states, some times lower continuum threshold values are 
favorable in order to completely eliminate the continuum above the resonance state. Therefore, in Fig. O 
we also include the result for y/so — 5.1 GeV. We see that we get a very narrow Borel window, and for 
values of the continuum threshold smaller than 5.1 GeV there is no allowed Borel window. Taking into 
account the variations on AP, sq and rric in the regions indicated above we get: 



mz^^_^ = (4.84 ±0.14) GeV , 



(20) 
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FIG. 5: The with — 1~ meson mass as a function of the sum rule parameter for different values of ^/sq: 
^/so = 5.1 GeV long-dashed line, ^/so = 5.2 GeV dashed line, ^/so = 5.3 GeV solid line and ^/so — 5.4 GeV 
dot-dashed line. The crosses indicate the region allowed for the sum rules. 



which is much bigger than the experimental value and bigger than the result obtained using the current 
with = Q- in Eq. (fH]) . 

For the value of the parameter A defined in Eq. ([5]) we get: 

A^^^_^ = (8.36 ±0.85) X lO^^ GeV^ (21) 

In conclusion, we have presented a QCDSR analysis of the two-point function of the recently observed 
Z+(4430) meson, considered as a tetraquark state, with a diquark-antidiquark configuration. Since the 
spin-parity quantum numbers of the Z^(4430) meson are not known, we have considered two different 
possibilities: — 0^ and = 1^. We have found a very goo d OPE convergence for these two cases, 
although this is not in general the case for tetraquark states [30] ■ We got a mass in some agreement 
with the experimental result in the case with — 0^. However, in the case = 1^, we got a much 
higher value for the mass. This is consistent with the expectation from the constituent quark model, since 
in this model the scalar diquark component in the color anti-triplet configuration is the most attractive 
quark-quark channel. 

Comparing our result, for the case — 0~, with the case where the Z+(4430) meson was considered 
by using a D*Di molecular current, also with = 0~ [§;], the differences are also not really big. Since 
there is no one to one corresponde between the structure of the current and the state, we can not conclude 
that the Z+(4430) is better explained as a molecular state than as a diquark-antidiquark state. However, 
comparing the results obtained for the quantum numbers = 0~ and 1^, from our calculations we 
conclude that the Z+(4430) is probably a — 0^ state. 
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